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Abstract. The observed properties of the close binary AE Aqr indicate that the mass transfer in this system operates via the 
Roche lobe overflow mechanism, but the material transferred from the normal companion is neither accreted onto the surface of 
the white dwarf nor stored in a disk around its magnetosphere. As previously shown, such a situation can be realized if the white 
dwarf operates as a propeller. At the same time, the efficiency of the propeller action by the white dwarf is insufficient to explain 
the rapid braking of the white dwarf, which implies that the spin-down power is in excess of the bolometric luminosity of the 
system. To avoid this problem we have simulated the mass-transfer process in AE Aqr assuming that the observed braking of 
the white dwarf is governed by a pulsar-like spin-down mechanism. We show that the expected Ha Doppler tomogram in this 
case resembles the tomogram observed from the system. We find that the agreement between the simulated and the observed 
tomograms is rather good provided the mean value of the mass-transfer rate < M >~ 5 x 10'* g s"' . Three spatially separated 
sources of Ho- emission can be distinguished within this approach. The structure of the tomogram depends on the relative 
contributions of these sources to the Ha emission and is expected to vary from night to night. 
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1. Introduction 

AE Aquarii is a peculiar nova-like star at a distance of ~ 100 + 
30 pc (Welsh et al. il 9951 Friediung |1997> . It is a non-eclipsing 
, close binary system with an orbital period Porb ~ 9.88 hr and 
• orbital eccentricity e « 0.02 (Chincarini & Walker 198 li. The 
normal companion (secondary) is a K3-K5 red dwarf on or 
close to the main sequence (Bruch .1991. Welsh et al. .1995 1. 
' The primary is a magnetized white dwarf rotating with the pe- 
riod ~ 33 s (Patterson 1979 Eracleous et al. [D^ . The 
incUnation angle of the system and the mass ratio are lim- 
ited to 50° < / < 70°, and 0.58 ^ (q = M2/M1) < 0.89, 
respectively, and the mass of the white dwarf is evaluated as 
0.8 ^ Ml < IMq (Reinsch & Beuermann [T994l Welsh et al. 
[T995t . 

The system emits detectable radiation in almost all parts of 
the spectrum. It is a powerful non-thermal flaring radio source 
(Bastian et al. 119881 Meintjes & Venter 2003 and references 
therein) and, possibly, the source of very high energy y-rays 
(Bowden et al. [T99 2, Meintjes et al. 119941 see, however, Lang 
et al. 199S). The optical, UV, and X-ray radiation of the system 
is predominantly thermal and comes from at least three dif- 
ferent sources. The visual light is dominated by the secondary 
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(Bruch lTWn Welsh et aL fD^ . The contribution of the pri- 
mary is observed mainly in the form of 33 s (and 16.5 s) co- 
herent oscillations detectable in the optical, UV, and X-rays 
(Patterson [DT^ Eracleous et al. [T^ Patterson et al. lTMJt . 
The remaining light comes from a highly variable extended 
source, which manifests itself in the blue/UV continuum, the 
optical/UV broad single-peaked emission lines, and the non- 
pulsing X-ray component. This source is associated with the 
mass-transfer process and is suspected of being responsible for 
the peculiar rapid flaring of the star (for discussion see e.g. 
Eracleous & Horne [T996t . 

AE Aqr is currently assigned to the DQ Her subclass of 
magnetic Cataclysmic Variables (CVs). The members of this 
subclass are interacting low-mass close binaries, in which the 
degenerate companions are magnetized white dwarfs rotating 
with periods <s: Porb and accreting material from a Keplerian 
disk (see e.g. Warner [T995> . However, extensive investigations 
during the last decade have clearly shown that AE Aqr does not 
fit in this model. Namely, studies of the 33 s pulsations in the 
optical/UV (Eracleous et al. ll994> and X-rays (Reinsch et al. 
IT^ Clayton & Osborne [T9951 Choi et al. revealed 
that the contribution of the white dwarf to the system radia- 
tion is significantly smaller than previously assumed within the 
accretion-powered white dwarf model. Furthermore, analysis 
of the Ha Doppler tomogram of AE Aqr has shown no evi- 



2 



N.R. Ikhsanov et al.: On the mass transfer in AE Aquarii 



dence of an accretion disk in the system (Wynn et al. 119971 
Welsh et al. 1998 1. Finally, de Jager et al. J1994> reported a 
mean spin-down rate of the white dwarf Pq - 5.64xl0"''^ss"'. 
The remarkable stability of the observed braking over the span 
of 14.5 yr suggests that the entire white dwarf is spinning down 
(for a detailed discussion see Welsh 1999). This allows us to 
evaluate the spin-down power of the white dwarf as 

L,d = 6 X 10^^ /5o P-^i (P/Po) erg s"', (1) 

where /50 and f 33 are the moment of inertia and the spin pe- 
riod of the white dwarf expressed in units of 10^"gcm^ and 
33 s, respectively. L^a exceeds the luminosity of the system ob- 
served in the UV and X-rays by a factor of 120-300 and even 
the bolometric luminosity by a factor of more than 5 (hereafter, 
the distance to the star is adopted as lOOpc). This indicates 
that the spin-down power dominates the system energy budget 
and raises a question about the form in which this energy is 
released. 

Among possible answers to this question the following two 
are currently under discussion. The first, a so called 'magnetic 
propeller' model, was presented by Wynn et al. (1997), who 
suggested that the rotation rate of the white dwarf deceler- 
ates by means of interaction between its fast rotating magneto- 
sphere and the material inflowing from the secondary. The sec- 
ond, a so called 'pulsar-like white dwarf model, was suggested 
by Ikhsanov (1998), who indicated that the observed braking 
of the white dwarf could be explained in terms of the canonical 
pulsar- like spin-down mechanism (Pacini 1968; Goldreich & 
Julian [R69), provided its surface magnetic field is 50 MG. In 
this paper we address the comparative analysis of these models. 
The basic statements of these approaches are briefly discussed 
in the following two sections. In Sect. 4 we present the results 
of our simulation of the Ha Doppler tomogram of AE Aqr. The 
adopted assumptions are summarized in Sect. 5. The basic con- 
clusions are given in Sect. 6. 

2. Propeller action by the white dwarf 

The observed properties of the optical/UV emission lines un- 
ambiguously indicate that a relatively intensive mass-transfer 
takes place between the system components of AE Aqr In par- 
ticular, the narrow component of the Balmer emission lines is 
nearly in anti-phase to the absorption lines of the red dwarf 
(Reinsch & Beuermann ll994t . that suggests its origin is near 
the white dwarf. Furthermore, the evaluated velocity and lumi- 
nosity of the emission line source are significantly larger than 
those typically expected in the wind of red dwarfs (see e.g. 
Eracleous et al. '1994'). These properties speak in favor of an 
association of the emission line source with the material trans- 
ferred from the normal component through the Roche lobe of 
the white dwarf. 

The rate of mass transfer in AE Aqr is still a subject of 
discussion. A lower limit to this parameter can be derived as- 
suming that the radiation of emission lines is generated inside 
the Roche lobe of the white dwarf and is powered by the accre- 
tion energy. In this case one finds M > 10'^L3i7?ioMQg g s"', 
where L31 is the luminosity of the emission line source (see e.g. 
Table 3 in Eracleous & Horne .1996.) . and M0.9 is the mass of the 



white dwarf in units of 0.9 Mq. Riq is the distance of the closest 
approach of the material responsible for the observed emission 
lines to the white dwarf expressed in units of 10'*^ cm. This pa- 
rameter can be limited using the expression R $ GM„i/V^^, 
where V^m is the velocity of the emitting material derived from 
the observed width of the emission lines. 

The above limit to M represents the minimum possible 
value of the mass-transfer rate in AE Aqr and, as will be shown 
below, is significantly underestimated. Nevertheless, this esti- 
mate plays an important role in the identification of the mass- 
transfer mechanism. Indeed, the derived value exceeds the 
maximum possible rate of mass capture by the white dwarf 
from the wind of its companion by more than three orders 
of magnitude (Ikhsanov 19971. This justifies that the mass- 
transfer in AE Aqr operates via the Roche lobe overflow mech- 
anism and hence, the secondary overflows its Roche lobe and 
loses material through the LI point towards the primary. 

However, a relatively low X-ray luminosity of AE Aqr 
(Lx ~ 10^' ergs"', see e.g. Choi et al. .1999 ) and the structure 
of the Ha Doppler tomogram derived by Wynn et al. ( I1997> 
and Welsh et al. (1998) argue against the possibility that the 
material transferred from the red dwarf is either accreted onto 
the surface of the white dwarf or stored in a disk. To solve this 
paradox the hypothesis has been invoked that the material flow- 
ing into the Roche lobe of the white dwarf is ejected from the 
system without forming a disk. 

An effort to reconstruct the mass-transfer picture within this 
hypothesis was first made by Wynn et al. Cl997j . They mod- 
elled the stream as a set of diamagnetic blobs, which move 
through the fast rotating magnetosphere of the primary, inter- 
acting with the local magnetic field via a surface drag term. 
In this case, the trajectories of the blobs differ from the bal- 
listic case, and the stream is able to leave the system with- 
out forming a disk if the magnetic moment of the primary 
is fj. lO-'^Gcm-'. In particular, putting jj ^ lO-'^Gcm'' 
(i.e. within the expected range of the magnetic moments of 
Intermediate Polars), Wynn et al. (1997 ) found that blobs reach 
the escape (maximum) velocity of Vesc ~ 1000 km s"' at the 
closest approach to the white dwarf, ro ^ 10'" cm, and leave 
the system without forming a disk with an average velocity 
Voo ~ 300km s '. The ejection of material in this scenario oc- 
curs due to propeller action by the white dwarf, which is also 
assumed to be responsible for the observed braking of the pri- 
mary. 

The Ha Doppler tomogram calculated within this model is 
similar to the tomogram observed in AE Aqr in several im- 
portant aspects. In particular, neither shows azimuthal sym- 
metry, and the emission is not centered on the white dwarf 
but is primarily in the lower-left quadrant (Vx, Vy both nega- 
tive). These similarities indicate that the picture reconstructed 
by Wynn et al. (1997 ) is qualitatively correct. 

At the same time, some of the quantitative predictions of 
the 'magnetic propeller' model have not been observation- 
ally confirmed. As shown by Welsh et al. (1998), the ob- 
served tomogram does not contain the high velocity 'loop' 
{\V\ ~ 700 - 1000 km s-') predicted by Wynn et al. ( fT997l 
see Fig. 3), and on the other hand it shows that the contribu- 
tion of material at low velocities {\V\ S 100 km s"') is signif- 
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icantly larger than expected from the simulated picture. These 
inconsistencies forced Welsh et al. J 19981 1 to suggest that the 
heating of the blobs, when they pass the acceleration region 
at their closest approach to the white dwarf, is negligible, and 
therefore that their contribution to the Ha emission of the sys- 
tem is small. Following this assumption they have placed the 
region of energy release outside the Roche lobe of the pri- 
mary where the trajectories of blobs of different masses cross 
each other and collisions of the ejected blobs can be expected. 
However, as mentioned by Welsh ( 1999). some of the prop- 
erties of AE Aqr (such as the large velocities in the emission 
lines during flares and the existence of high-excitation emission 
lines) remain puzzling even in this so called 'colliding blobs' 
scenario. 

Another difficulty with the 'magnetic propeller' model 
has been mentioned by Ikhsanov U998J . Meintjes & de Jager 
(EOOO), and Choi & Yi SKMi . As fliey have shown, the efli- 
ciency of the propeller action by the white dwarf under the con- 
ditions of interest is not sufficient to explain the observed rapid 
braking of the primary. Indeed, following Wynn et al. ( I1997t 
one could assume that almost all spin-down power of the white 
dwarf is transferred into the kinetic energy of the ejected ma- 
terial. However, the kinetic luminosity of the ejected blobs is 
obviously limited to 

Lkin ^ {l/2)MVi, - 5 X 10^2^17 Vi erg s-\ (2) 

where = Vesc/lO^^cms"' and Mn is the mass transfer rate 
expressed in units of 10'^ gs '. Hence, for this assumption 
to be satisfied, the mass transfer rate in the system should be 
in excess of 10'^ gs"', which is inconsistent with the value 
of M derived from observations (see e.g. Eracleous & Horne 
I1996t . On the other hand, the assumption that a significant part 
of the spin-down power is transferred into the thermal energy 
of the ejected gas contradicts the derived value of the ratios 
Luv/i'sd ~ Lji/L^d ^ 1- Therefore, the question about the 
nature of the observed braking of the white dwarf within the 
'magnetic propeller' model remains open. 

The problems mentioned above indicate that the 'magnetic 
propeller' model cannot provide us with a complete picture of 
AE Aqr, and an improvement of this model is required. They 
also show that the problem that should be addressed first in 
any further improvements is the spin-down mechanism of the 
white dwarf. As long as this problem remains unsolved, the 
form in which the spin-down power is released turns out to be 
unclear, and therefore, the major part of the energy released in 
the system is not taken into account. 

At the same time, for a solution of the spin-down prob- 
lem to be reliable it should also meet the diskless mass-transfer 
criteria. In this light, the improvement suggested by Meintjes 
& de Jager (^000 1 cannot be applied to the interpretation of 
AE Aqr, since their approach requires the existence of a clumpy 
disk around the white dwarf. On the other hand, the model of 
Choi & Yi (2000), in which the spin-down power is assumed 
to be spent in the emission of gravitational waves, cannot be 
accepted either Although the mass-transfer picture within this 
model is similar to that reconstructed by Wynn et al. J1997> . the 
adopted mass distribution over the primary surface is unreliable 
(for a detailed discussion see Ikhsanov & Beskrovnava .2002< . 



Among the improvements of the 'magnetic propeller' model so 
far discussed in the literature, only the pulsar-like spin-down 
model meets the criteria. The reliability of this improvement is 
discussed in the following section. 

3. Pulsar-like spin-down 

The hypothesis of pulsar-like spin-down of the white dwarf 
in AE Aqr has a certain observational basis. A situation in 
which the spin-down power of a star exceeds its bolometric 
luminosity significantly is unique for CVs as well as for all 
presently known accretion-powered sources. At the same time, 
this situation is typical for the spin-powered pulsars, whose lu- 
minosity constitutes only a small fraction (~ lO^^* - 10"') of 
their spin-down power (for a review see Manchester & Taylor 
119771 Hartmann 1 19951 1. Furthermore, while the appearance of 
AE Aqr in X-rays is very atypical for the accretion-powered 
white dwarfs (Clayton & Osborne '1995), it resembles the ap- 
pearance of spin-powered pulsars observed in the ROSAT en- 
ergy range (see e.g. Becker & Triimper 1 19971 1. For instance, 
the X-ray spectrum is significantly softer than those typically 
observed from accretion-powered compact stars, and the ratio 
of the luminosity of the pulsing component to the spin-down 
power is close to lO^^*. Finally, as reported by Meintjes et al. 
( fl994) , the intensity of the very high energy y-ray emission 
detected from AE Aqr significantly exceeds the intensity of ra- 
diation emitted in other parts of the spectrum. Such behavior is 
also typical for spin-powered pulsars (see e.g. Tompson 1996.I 
and is consistent with modern views on the processes of energy 
release in these sources (for a review see Michel' 1991). Thus, 
the investigation of the hypothesis that the braking of both the 
white dwarf in AE Aqr and the spin-powered pulsars is gov- 
erned by the same mechanism appears to be quite reasonable. 

As shown by Ikhsanov ( 1998 ), for this hypothesis to be ef- 
fective the dipole magnetic moment of the white dwarf should 
be as large as 

/i- 1.4x10^^2 ( ^ X^Gcm\ (3) 

^ ^H6x 1033 ergs- V ^ 

This implies that the mean strength of the magnetic field at 
the surface of the white dwarf is B{R^±) ^ 50 MG. Under 
these conditions the spin-down power is spent in the gener- 
ation of magneto-dipole waves and particle acceleration and 
hence is released mainly in undetectable parts of the spectrum. 
Therefore, the observed inequality Lboi < isd turns out to be 
naturally explained within this approach. 

The limitation of the magnetic field to 50 MG is consis- 
tent with present views on possible values of the surface field 
strength of white dwarfs (see e.g. Jordan 2001 ) . In particular, 
the magnetic field of white dwarfs in Polars is of the same order 
of magnitude (Cropper 1990 ; Chanmugam 1992'). However, it 
is significantly above the previous limit to the strength of the 
magnetic field of the white dwarf in AE Aqr derived by Bastian 
et al. ( 198? ) and Stockman et al. (1 9921 1 from the analysis of the 
circularly polarized optical emission. 

The reason for this inconsistency has recently been investi- 
gated by Ikhsanov et al. (2002). As they have shown, the lim- 
itation presented by Bastian et al. (|1988^ and Stockman et al. 
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<199 2t is model-dependent and is based on the assumption that 
the radiation of the white dwarf is powered mainly by the accre- 
tion of material onto its surface. However, the investigations of 
AE Aqr in the UV (Eracleous et al l 19941 1 and X-rays (Clayton 
& Osborne [T995l Choi et al. 1995Jl have clearly shown that this 
assumption is not valid. As presently recognized, the contri- 
bution of the hot polar caps to the visual radiation of AE Aqr 
does not exceed 0. 1 %-0.2% . In this situation the hot polar caps 
cannot be the source responsible for the circularly polarized ra- 
diation observed from the system. Otherwise, the intrinsic po- 
larization of the source proves to be in excess of 100%, that is 
obviously impossible (see Ikhsanov et al. 2002 ). Therefore, the 
above mentioned inconsistency cannot be used as an argument 
to reject the possibility of the white dwarf in AE Aqr having a 
magnetic field as strong as 50 MG. 

A possible history of AE Aqr within the pulsar-like model 
is more complicated than that usually modelled within the 
'magnetic propeller' approach (Meintjes 2002 Schenker et al. 
12002 1. Indeed, a white dwarf with mass 0.9 M© and surface 
magnetic field 50 MG can only be spun up to the period of 33 s 
if the mass transfer rate during a previous epoch was in ex- 
cess of the Eddington limit (~ 3 x 10^' gs"') by a factor of 
3. Accretion with these characteristics resembles the process 
of the merging of a white dwarf with another star rather than 
the mass exchange between a main sequence red dwarf and a 
white dwarf in a close binary. The formation of a fast rotat- 
ing, strongly magnetized white dwarf due to its merging with a 
companion has been already discussed by Paczyriski (.1990 1. 
Following this scenario, one should assume that the white 
dwarf in AE Aqr is a product of the merging of a magnetized 
white dwarf and, possibly, a brown dwarf of mass ~ O.O3M0. 

However, as pointed out by Ikhsanov ( I1999t . the process 
of merging is not the only possible solution. An alternative ex- 
planation is based on the scenario of magnetic field amplifica- 
tion in very fast rotating compact stars ( Kluzniak & Ruderman 
11998 ' Spruit T999 1. According to Chanmugam et al. ( 1987> . 
the rotation of a white dwarf with mass Mwd = 0.9 Mq be- 
comes significantly non-uniform as its period decreases below 
Per - 20 s. The magnetic field inside the star in this state is 
winding up to ~ 10^ G on a time scale of a month, and mani- 
fests itself at the surface due to the buoyant instability produc- 
ing a surface field of ~ 10^ G. This allows us to envisage a situ- 
ation in which the magnetic field of the white dwarf in AE Aqr 
was amplified to its present value during the last stage of a pre- 
vious accretion-driven spin-up epoch. As shown by Ikhsanov 
J1999J . for this scenario to be eff'ective one has to assume that 
the initial magnetic moment of the white dwarf was 

A^o S 2 X 10^' Mjf Pll^' m'J^ gs-\ 

which implies that in a previous epoch AE Aqr was an ordinary 
member of the DQ Her subclass of CVs. 

Although both of the above mentioned scenarios lead to a 
rather complicated history of AE Aqr, it is clear that a solution 
of this problem within the pulsar-like spin-down model is not 
impossible. A more precise investigation, however, is not effec- 
tive as long as the presently observed stage of the system is not 
well identified. Therefore, in this paper we will focus mainly 
on the analysis of currently observed properties of AE Aqr. 



The natural solution of the spin-down problem is not the 
only advantage of the pulsar-like model. It also gives a rea- 
sonable explanation of some properties of AE Aqr observed 
in the high-energy parts of the spectrum. In particular, it pre- 
dicts the maximum energy of particles accelerated by the white 
dwarf fi™" ~ 2 X 10'^ eV (see Eq. 10 in Ikhsanov fl^ . This 
prediction is consistent with the characteristics of the TeV y- 
ray events observed from the system (Meintjes et al. I1994> . 
Furthermore, it also allows us to associate the origin of the puls- 
ing UV and X-ray emission with the processes of non-thermal 
energy release in the magnetosphere of the white dwarf, i.e. 
particle acceleration in the inner (and, possibly, outer) gap and, 
correspondingly, the impact of particles responsible for the 
back-flowing current onto the surface of the white dwarf at the 
magnetic pole regions (for a discussion see Ikhsanov 1998 and 
references therein). The latter process should lead to the heat- 
ing of the surface of the white dwarf, while the radiative losses 
of relativistic electrons are expected to be observed in the X-ray 
part of the spectrum. Under these conditions, the luminosity of 
pulsing emission in the UV would be comparable with that of 
the pulsing X-ray component, and the area of the hot polar caps 
can be limited to (see Eq. 11 in Ikhsanov . 1998.) 

Ape < 3 X 10'" Tj,\, Ml'l^ n'l^^ Rl, Mo.9 cm^, (4) 

where //34.2 and T^s.s are the magnetic moment and the radius of 
the white dwarf expressed in units of 10^"*^ G cm^ and 10** cm, 
respectively. 770.37 - 77/0. 37 is the parameter accounting for the 
geometry of the accretion flow, which in the case of a stream is 
normalized following Hameury et al. (I1986> . 

Both of these predictions are in good agreement with corre- 
sponding properties of AE Aqr (see e.g. Eracleous et al. 19941 
Choi et al. ll999> . and they allow us to avoid a very controver- 
sial assumption about the accretion nature of the pulsing UV 
and X-ray radiation of AE Aqr (for a discussion see Choi et al. 
[T999l and Ikhsanov l200il . 

Finally, the pulsar-like model naturally leads to a conclu- 
sion about the diskless mass transfer in the system. Indeed, 
within this approach the Alfven radius of the white dwarf, 

Pa - 3 X 10'" 770.37 1^122 ^n" '^^i' (5) 

is larger than the circularization radius for all reasonable values 
of M. This means that diskless mass transfer in the system is re- 
alized, and moreover, this is expected independently of whether 
the stream transferred through the LI point is initially homoge- 
neous or inhomogeneous. In this case the assumption that the 
stream disintegrates into blobs before the point of its closest 
approach to the white dwarf (as adopted by Wynn et al. ll997> 
turns out to be unnecessary. Instead, one can envisage a sce- 
nario in which the stream disintegration occurs in the vicin- 
ity of the Alfven surface of the white dwarf, where such a 
disintegration, according to Arons & Lea (1980), is expected. 
However, is the Hff Doppler tomogram expected within this 
scenario consistent with the tomogram observed from the sys- 
tem ? The analysis of this question is addressed in the next sec- 
tion. 
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4. Simulation of Ha Doppler tomogram 

We consider a situation in which the secondary star overflows 
its Roche lobe and loses material in the form of a stream 
through the LI point. The stream flows into the Roche lobe 
of the white dwarf at the speed of sound and initially follows 



a ballistic trajectory. Following Wynn et al. J1997> . we assume 
that at a certain point (its location depends on the scenario con- 
sidered) the stream disintegrates into a set of large diamagnetic 
blobs. The blobs interact with the magnetospheric field of the 
white dwarf via the drag term. Due to this interaction their tra- 
jectories are modified by the magnetic acceleration 



§mag = -k[V - Vflx, 



(6) 



where k ~ B^/cAPbk is the drag coefficient, V and Vf are the 
blob and field velocities, and the suffix ± denotes the velocity 
component perpendicular to the field lines, ca is the Alfven 
speed in the interblob plasma, which under the conditions of 
interest can be approximated by the speed of light, pb and l\, 
are the density and radius of the blobs. 

As shown by Wynn & King \l995j, the drag coefficient can 
be expressed in the form k ~ koir/ro)'", where ko, n, and ro 
are constants. Setting phiro) - 10""p_ii gem""* and - 
lO^/gcm we evaluate the parameter fco as (for discussion see 
Wynnetal. 1997 1 



^0-3.3x10-'' B\ra)p-_\, Ig' s" 



(7) 



Simulating the trajectory of the material within the Roche 
lobe of the white dwarf we have assumed that the drag interac- 
tion between the stream and the magnetic field of the primary 
before the point of the stream disintegration is small and there- 
fore, the stream before this point follows a ballistic trajectory. 
This assumption is reasonable if the point of the stream disinte- 
gration is located at a distance r ^ Rdm, where Rare is the circu- 
larization radius, which is about 2.5 X 10'° cm for the parame- 
ters of AE Aqr (see Eq. 4. 17 in Frank et al.'1985V Starting with 
the point at which the stream disintegrates into blobs the drag 
term becomes important and therefore, the trajectories of the 
blobs deviate significantly from the ballistic ones. Following 
this method, we have performed several runs of calculations 
placing the point of disintegration at the LI point, at the point 
of the closest approach of the stream to the white dwarf, ro, and 
at a few intermediate points located within the interval [LI, ro]. 
This scheme of calculations has been applied to both the 'mag- 
netic propeller' and the 'pulsar-like white dwarf approaches. 
The results of our simulations are presented in the following 
two subsections. 

4.1. The magnetic propeller approach 

To test the consistency of the results of our calculations with 
those previously derived by other authors we have simulated 
the trajectories of blobs using the parameters of AE Aqr 
adopted by Wynn et al. (1997") as follows. System parameters: 
mass ratio q - 0.64, orbital period Poib - 9.9 hr, and inclina- 

o 

tion angle / = 55 , and parameters of the white dwarf: mass 
Ml - Q.9Mq, dipole magnetic moment jj = lO^'^Gcm^, and 
spin period f s = 33 s. The stream of material transferred from 



the secondary has been modelled as a set of diamagnetic blobs. 
The distance to the point of the stream disintegration has been 
assumed to satisfy the condition rdis » Rchc- The radius and 
the density of blobs at their closest approach to the white dwarf 
have been taken as Zb = 10^ Ig cm and pb - 10"" p-n g cm""*, 
respectively. The value of the parameter n has been chosen to be 
n - 2. Finally, the parameter ro has been limited to ro t R^ia, 
where 



™"" lo.64j ll.8xlO>icmj 



cm 



(8) 



is the distance of the closest approach of a homogeneous stream 
to the white dwarf (see Eq. 2.14 of Warner il995< ). 

The trajectories of the blobs and the expected Ha Doppler 
tomogram simulated under these conditions are shown in 
Fig-E The best fit to the observed tomogram is found for 
ko ^ (0.8 - 1.3) x 10"^ s"'. Lines 1^ represent the ti-ajecto- 
ries of blobs of different mass with the mass of the blobs de- 
creasing from line 1 to line 4. The more massive the blob the 
smaller the distance to which it approaches the white dwarf. 
The trajectories of blobs of different masses intersect beyond 
the Roche lobe of the white dwarf and, therefore, collisions of 
the ejected blobs in this region can be expected. The location 
of the collision zone is shown in panels a and b of Fig.^as a 
hatched region. This region represents the expected structure of 
the Ha Doppler tomogram within the 'colliding blobs' model. 
The structure of the tomogram derived within the approach of 
Wynn et al. ( 119971 . i.e. under the assumption that the Ha emis- 
sion comes mainly from blobs moving through the Roche lobe 
of the white dwarf, is shown in panel c. 

As is easy to see, the derived tomograms are almost identi- 
cal to those presented by Wynn et al. ( 1997 see Fig. 3 of their 
paper) and Welsh et al. ( 1998 see Fig. 14 of their paper). This 
proves that our code is working properly and can be used for 
further analysis. 

4.2. The pulsar-like white dwarf approach 

The simulation of the stream trajectory within the 'pulsar- like ' 
model differs from that in the frame of the 'magnetic propeller' 
model in several important aspects. First, the dipole magnetic 
moment of the white dwarf within the 'pulsar- like ' model is 
assumed to be yu ^ 1.4 x lO^'^Gcm^, i.e. a factor of 100 larger 
than that adopted within the 'magnetic propeller' model. 
Second, the value of the parameter ro is limited to 



>Ra^3x 10'° 770.37 ^u'' ^0.9^' '^m- 



(9) 



This limitation reflects the fact that the magnetic field pressure 
at the Alfven radius is strong enough to prevent the stream from 
approaching the white dwarf closer than R\. Within the 'mag- 
netic propeller' model Ra < Rmin, and hence the pressure by 
the dipole field of the primary in the radial direction (within 
the frame centered at the white dwarf) can be neglected. But 
in the case of the 'pulsar-like' model the Alfven radius of the 
white dwarf is larger than Rmm- This means that the radial ve- 
locity of the inflowing material at Ra rapidly drops to zero with 
corresponding heating and possibly disintegration of the initial 
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Fig. 1. Trajectories of blobs in coordinate a and the velocity b and c space within the "magnetic propeller" approach. The position 
of the white dwarf is marked by the cross. The Roche lobe is shown for both system components in panel a, and for the normal 
component in panels b and c. The system separation is used as a unit scale in panel a. The lines 1-4 represent the trajectories of 
blobs of different masses in descending order with density. The hatched region in panels a and b indicates the location of the 
collisions zone of the blobs (see text for further details). 



stream. The further trajectory of the material depends on the ef- 
ficiency of the drag interaction between the blobs and the mag- 
netospheric field of the primary. If this interaction is effective 
enough for the blobs to be accelerated to the escape velocity in 
the azimuthal direction, they will leave the system. Otherwise, 
the material will follow an almost circular orbit around the 
magnetosphere of the white dwarf. 

Finally, the assumption about the stream disintegration at 
the LI point within the 'pulsar-like' model is not necessary. 
The point at which the stream disintegrates into diamagnetic 
blobs in this case can be located anywhere between the LI point 
and the Alfven surface of the primary. As mentioned above, in 
both cases (homogeneous and inhomogeneous stream) diskless 
mass transfer is expected within this model. Therefore, paral- 
lel to the traditionally considered case of an inhomogeneous 



stream at the LI point we also simulated the mass transfer as- 
suming that the point of disintegration of the initially homoge- 
neous stream is located at the Alfven surface of the primary. 

Evaluating the structure of the Ha Doppler tomogram, one 
should also take into account that a third source of Ha emis- 
sion (in addition to the stream passing through the magneto- 
sphere of the white dwarf and the region of blob collision) 
can be expected within the pulsar-like model of AE Aqr This 
source is associated with the region where the magneto-dipole 
waves emitted by the white dwarf are absorbed by the back- 
ground material surrounding the system. According to Rees & 
Gunn ( 1974) , the distance to this region, /?abs, can be found by 
equating the pressure of the magneto-dipole radiation, p^^ - 
LsdlA-ncR^, with the thermal pressure of the surrounding mate- 
rial, ppi = (l/2)pooVi^. The interaction between the waves and 
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Fig. 2. Trajectories of blobs in the coordinate {upper panels a-d) and the velocity {bottom panel) space within the pulsar-like 
model of AE Aqr. The position of the white dwarf is marked with a black dot. Panels a, b, and c were calculated for mass transfer 
rates of 10'^ g s"', 5 X 10'^ g s and lO'^ g s ' , respectively. The Unes in these panels show the trajectories of blobs of different 
masses. Panel d is the superposition of panels a, b, and c taken with similar weights. The bottom panel represents the expected 
structure of the Ha Doppler tomogram of the system within this approach. The calculations were made on the assumption that 
the stream is disintegrated into the blobs at the LI point. For a detailed description see text. 
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Fig. 3. The same as Fig|2 except for the assumption that the point of the stream disintegration is located at the Alfven surface of 
the white dwarf. 
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the gas leads to the formation of a shock in which the energy of 
waves is converted into the thermal energy of plasma, radiation, 
and accelerated particles. 

To estimate /?ab,s in AE Aqr, we have taken into account 
that the circumbinary medium of this system is contributed to 
mainly by the material ejected due to the propeller action of 
the white dwarf. As shown by Wynn et al. (tl99 7 1, this material 
flows out within the orbital plane of the system following a spi- 
ral. The position angle of the spiral, however, changes with the 
orbital motion of the system. Therefore, the distribution of the 
material surrounding the system has an azimuthal symmetry. 
Simulation of the stream-like ejection in the rotating system has 
shown that a circumbinary disk-like envelope with inner radius 
/?cbe S 5 X 10" cm forms around the system. The mean velocity 
of the ejected material at this distance is yb(^cbe) ~ 100 km s"', 
and therefore the derived scale is comparable to PoibVb(/?cbe)- 

The thickness of the envelope is determined by the thermal 
expansion of blobs, and its mean value can be normalized as 
•Zo(^cbe) - 10"^ Zio cm. This allows us to limit the mean density 
of the envelope material to 



At s 10-'" ll Z 



14 ;3 7-3 
10 



Phiro) 
10-"gcm- 



gcm" 



(10) 



Under these conditions the distance at which the magneto- 
dipole waves are absorbed by the material of the envelope can 
be evaluated as 



TJabs ^ 2 X lO'VllfV,' 



6 X 10^^ erg s" 



1/2 



cm. 



(11) 



At this distance the circumbinary envelope occupies about 
0.5% of the area of a sphere, and the velocity of the material 
lies within the interval ~ 50 - 150 km s"'. This indicates that 
the envelope will contribute to the Doppler tomogram of the 
system at low velocities with intensity $ 5 x lO^^^Lsd- 

The results of the simulation of the stream trajectories 
within the pulsar-like model of AE Aqr are shown in Figs.|2] 
and|3 These figures differ in the basic assumption about the lo- 
cation of the region of the stream disintegration to the diamag- 
netic blobs. Namely, in the first run of the calculations (Fig.|5J 
this region has been placed at the Ll point, and in the second 
run we have assumed that the stream disintegration occurs at 
the Alfven surface of the primary (Fig.|3}- The values of the 
system parameters (except for ji and ro) in all calculations were 
chosen to be the same as those adopted by Wynn et al. ('1997i. 

In both runs the calculations were made for three differ- 
ent values of the mass transfer rate: M - lO'^gs ' (panel a), 
M = 5 X 10'^ gs"' (panelb), and M = 10'^ gs"' (panelc). The 
value of the parameter rf was chosen as 0.5. Panel (d) shows 
the mean picture of the mass transfer on the time scale of the 
orbital period of AE Aqr. This picture is derived by taking a 
superposition of states (a), (b), and (c) with equal weights. The 
expected mean structure of the Hor Doppler tomogram derived 
in the first and the second runs are presented at the bottom of 
Fig.|2]and Fig.|3j respectively. 

One finds the best agreement between the simulated and 
the observed tomogram for < M >- 5 x 10'* gs"', and for 
ko — 4.8 X 10"^ s"' (assuming the stream to be disintegrated 



at Ll) and - 1 A x 10"^ s"' (assuming the stream to be 
disintegrated at R\). Here < M > denotes the average mass- 
transfer rate on the time scale of the orbital period. Under these 
conditions the average value of the Alfven radius of the white 
dwarf is R\ ^ 4.9 x 10'°cm. This allows us to evaluate the 
average bolometric luminosity of the stream passing through 
the magnetosphere as Lst ~ (< M >)GM^d/RA - 10^^ ergs"'. 
This value slightly exceeds the luminosity of the low-velocity 
source situated beyond the light cylinder of the white dwarf 
and is close to the average luminosity of the extended com- 
ponent of radiation in AE Aqr evaluated by van Paradijs et al. 
( 1989 1, and Eracleous & Horne ( I1996t from the optical and UV 
observations. The main features of the derived tomograms are 
discussed in the following section. 



4.3. Comparative analysis of ttie derived tomograms 

The tomograms calculated within the 'pulsar-like' and the 
'magnetic propeller' models have a number of similarities. In 
particular, neither shows azimuthal symmetry, and the emission 
is not centered on the white dwarf but is primarily in the lower- 
left quadrant (V^, Vy both negative). Because of these proper- 
ties all of the simulated tomograms resemble the Doppler to- 
mogram observed from AE Aqr. 

There are, however, several important differences. First, 
the upper limit to the velocity of the stream at the clos- 
est approach to the white dwarf within the pulsar-like model 
is smaller by a factor of 2 than that within the 'magnetic 
propeller' approach. This means that the emission associated 
with the stream passing through the magnetosphere within the 
pulsar-like model is produced at velocities < 500km s '. Since 
the blobs have different masses, dispersion of their veloci- 
ties occurs. Furthermore, the velocity of blobs at their clos- 
est approach to the primary depends on the mass-transfer rate. 
Superposition of these effects leads to a situation in which the 
contribution of the stream at ro appears in the mean Doppler 
tomogram in the form of a spread loop, which is centered 
at (-250kms"'; -350kms"') and has the size of lAVxl ~ 
600kms-' and lAVyl ~400kms-'. 

The emission at these velocities is present in the observed 
Ha Doppler tomogram of AE Aqr. This indicates that blobs 
passing through the magnetosphere within the pulsar-like ap- 
proach are expected to be hot, and their contribution to the Ha 
emission of the system is significant for all reliable values of 
M. Therefore, the problem of the 'missing radiation' from the 
loop associated with the blob trajectories (mentioned by Welsh 
et al. 19981 does not occur, as one assumes the surface mag- 
netic field of the white dwarf to be of the order of 50 MG. 

The second feature of the tomogram simulated within the 
pulsar-like approach is the significant contribution of the ma- 
terial situated beyond the light cylinder of the white dwarf. 
The radiation of this source is emitted at velocities ~ 50 - 
150 km s"' and is powered by the energy of the magneto-dipole 
waves. The luminosity of this source is almost independent 
of the variations of the mass transfer rate and is comparable 
with the luminosity of the stream passing through the mag- 
netosphere at M ~ 5 X lO'^gs '. The contribution of this 
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low-velocity source to the Ha system radiation is seen in the 
center of the Doppler tomogram as a spread spot of radius 
~ 150kms-i. 

An additional, intermediate-velocity source of Ha emis- 
sion is located at a distance of about 1-3 times the binary 
separation. This source is associated with the region of possi- 
ble blob collisions. Our simulations indicate that in both mod- 
els collision of blobs can occur and that the relative velocity 
of the colliding blobs is of the order of Vh-h ~ 100km s '. 
Assuming that all blobs are involved in the collision pro- 
cess, one can limit the rate of energy release in this region to 
Lb-b ^ 5 X 10^"Mi7(yb-b/100kms-')2ergs-i. This indicates 
that the contribution of this source within the pulsar-like model 
can be significant at relatively high mass-transfer rates but can 
hardly be recognized when the mass transfer rate drops below 
lO'^gs-i. 

Finally, our simulations have shown that the structure of the 
tomogram within the pulsar-like model is sensitive to variations 
in the average mass-transfer rate in the system. As < M > de- 
creases, the Alfven radius of the white dwarf becomes larger. 
In this case the material moving through the magnetosphere 
turns out to be ejected at smaller velocities, and its contribu- 
tion to the Ha emission of the system decreases. Therefore, 
the tomogram, under these conditions, is dominated by the 
low velocity component. If, however, < M > is large during 
the period of observations, the tomogram is dominated by the 
'high-velocity spot', which in this case appears in the lower-left 
quadrant at velocities 350 - 500 km s '. Hence, the observed 
night-to-night variations of the tomogram (see Fig. 10 of Welsh 
et al. 1199 8) can be interpreted within the pulsar-like model in 
terms of the variations of the mass-transfer rate in the system. 
The range of these variations implies changes in the efficiency 
of the propeller action by the white dwarf within the interval 
Lit/Lsd ~ 0.01 - 0.4, and therefore, its contribution to the ob- 
served braking of the white dwarf remains small. 

5. Discussion 

It is widely believed that the 'magnetic propeller' is the only 
approach which provides a plausible interpretation of the Ha 
Doppler tomogram observed in AE Aqr. Following this no- 
tion, almost all manifestations of the system during the past 
5 years have been discussed solely around the hypothesis that 
the spin-down of the white dwarf is governed by the propeller 
spin-down mechanism. 

However, as shown in this paper, the Ha Doppler tomo- 
gram expected within the pulsar-like approach also resembles 
the observed tomogram. Furthermore, the agreement between 
the expected and the observed tomograms within this approach 
turns out to be even better than within the 'magnetic propeller' 
model. As mentioned in Sect. 3, the basic assumptions adopted 
within the pulsar-like approach do not contradict any of the cur- 
rently observed properties of the system, but they allow us to 
invoke the models developed with respect to the spin-powered 
pulsars for the interpretation of properties which AE Aqr shares 
with at least several objects of this class. Therefore, an anal- 
ysis of the observed system properties within the pulsar-like 
approach appears to be quite reasonable. 



The present state of development of both the 'magnetic pro- 
peller' and the pulsar-like models is insufficient for recogniz- 
ing which of these approaches is more promising. The analy- 
sis of this question is beyond the scope of the present paper. 
Nevertheless, to clarify the basic statements of these models 
we summarize the assumptions adopted within each of these 
approaches. 

5.1. Assumptions adopted within the magnetic 
propeller approach 

The following 6 basic assumptions, currently adopted within 
the magnetic propeller model, can be distinguished: 

I. The dipole magnetic moment of the white dwarf is as- 
sumed to be of the order of lO-*^ G cm-^, i.e. within the expected 
range of the magnetic moments of Intermediate Polars. 

II. The spin-down power is assumed to be transferred pre- 
dominantly to the kinetic energy of the material ejected from 
the system due to propeller action by the white dwarf. As 
shown in Sect. 2.1.2, this assumption implies the mass-transfer 
rate in the system to be 

M>\0'^(—^]( ^ -Us-'. (12) 

ll0"'cm/\6x 103%-gs-V ^ ^ ^ 

III. The stream of material transferred from the secondary is 
assumed to be strongly inhomogeneous. Actually, this assump- 
tion implies that at least 99.99% of the material is transferred 
in blobs. Indeed, if the mass transfer rate of the homogeneous 
component exceeds 

/ ^ \-V/2 

Mbs ~ 10'^ ^ij^ M-f — g s-', (13) 

\"circ/ 

the homogeneous flow turns out to be able to reach the circu- 
larization radius and to form a disk around the magnetosphere 
of the white dwarf. 

It should be noted that the reason for such a strong inho- 
mogeneity is rather unclear. It might be connected with the 
magneto-flaring activity of the normal component or the beam 
instability in the region of the LI point. At the same time, it 
is unlikely that it can be explained in terms of the interaction 
between the stream and the magnetic field of the white dwarf 
since at distances /? » 7?a ffie energy density of the magnetic 
field is a few orders of magnitude smaller than ffie thermal en- 
ergy of the stream material. In particular, the solution of Axons 
& Lea ( 1980 ) derived for the regions R ~ R^is obviously not 
applicable in this case. 

IV. The temperature of blobs passing through the magne- 
tosphere is assumed to be small enough for their contribu- 
tion to the Ha emission of the system to be negligibly small. 
Otherwise, the presence of a high velocity loop associated with 
the trajectories of blobs interacting with the magnetic field of 
the white dwarf is expected (see Fig. 3 in Wynn et al. I1997> . 
Such a loop, however, is not seen in the observed tomogram 
(for discussion see Welsh et al. [T998l Welsh .1999) . 

V. Almost all the blobs expelled by the white dwarf are in- 
volved in a collision process as their trajectories intersect be- 
yond the system. Otherwise, the amount of hot material would 
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be insufficient to explain the observed Ho- emission. It should 
also be noted that for the observed flaring in the system to be 
associated with the collision of blobs, the mass transfer rate 
should be in excess of (see Sect. 4.3) 

where L33 is the luminosity of the flaring component ex- 
pressed in units of 10^^ ergs"' (see van Paradijs et al. 1989 
Beskrovnaya et al. 1996), and the value of Vb-b is normalized to 
the rising velocity of flares recently derived by Skidmore et al. 
J2003> from high-time-resolution spectroscopy of the system. 
Otherwise, the energy release in the region of blobs collision 
turns out to be smaller than the luminosity of strong flares ob- 
served from the system. 

VI. One has to assume that the blobs cool down very slowly 
(on a time scale of a few hours) or, for some reason, are heated 
again at larger distances from the system. Otherwise, the ori- 
gin of the low-velocity (^ 100 km s ') component of the sys- 
tem Ha emission, which is seen on the observed tomogram, 
becomes rather unclear (see e.g. Welsh 119991 Pearson et al. 
12003^. 

Within these assumptions a good agreement between the 
observed Ha Doppler tomogram and the tomogram simulated 
within the 'magnetic propeller' model can be achieved. 

5.2. Assumptions adopted within the pulsar-like white 
dwarf approach 

The basic assumptions of the pulsar-like model of AE Aqr are 
as follows: 

I. The dipole magnetic moment of the white dwarf is as- 
sumed to be as high as ju ^ 1.4 x lO^'^Gcm^. As shown in 
Sect. 3, this assumption implies a rather complicated scenario 
for the system evolution, which invokes processes poorly in- 
vestigated so far 

II. The stream of material transferred from the secondary 
is assumed to disintegrate into large diamagnetic blobs as it 
impacts onto the Alfven surface of the primary. 

III. It is assumed that the mass-transfer rate in the system 
varies by a factor of a few on the time scale of the orbital period. 
Due to these variations the structure of the tomogram changes 
from night to night. 

IV. We have also assumed that the stream does not strongly 
interact with relativistic particles accelerated in the potential 
gap situated in the vicinity of the white dwarf surface. This im- 
plies that the opening angle of the beam of particles accelerated 

o 

in the magnetic pole regions is ^ 18 (here, the value of is 
adopted as ^ 77 , see Eracleous et al. ll994> . If, nevertheless, 
there is interaction between the stream and the accelerated par- 
ticles, the energy balance in the material of the stream should 
be re-calculated. 

At the same time, none of the assumptions adopted within 
the 'magnetic propeller' approach is necessary in the frame 
of the pulsar-like model of AE Aqr. In particular, there is no 
strong lower limit to the mass-transfer rate, diskless mass trans- 
fer is expected independently of whether the stream inflowing 



through the LI point is homogeneous or inhomogeneous, the 
blobs passing the accelerating region are expected to be hot, 
and the low-velocity source is associated with the interaction 
between the magneto-dipole waves and the plasma expelled by 
the propeller action of the white dwarf. 

6. Conclusions 

We have shown that the Ha Doppler tomogram simulated 
within the pulsar-like white dwarf model of AE Aqr is sim- 
ilar to the observed tomogram in several important aspects. 
Namely, the emission is not centered on the white dwarf, it 
does not show azimuthal symmetry, and the strongest emis- 
sion occurs primarily in the lower-left quadrant at velocities 
< SOOkms-i. 

At least three sources of the Ha emission can be distin- 
guished within the considered approach: (1) the emission of the 
stream passing through the magnetosphere of the white dwarf 
(the high-velocity component: 350 - 500 km s"'), (2) the re- 
gion of blob collisions (the intermediate velocity component: 
200 - 300 kms"'), and (3) the region of interaction between 
the stream and the magneto-dipole radiation of the white dwarf 
(the low-velocity component: ^ 150 km s '). The relative con- 
tributions of these components to the system emission depend 
on the mass-transfer rate, and therefore the structure of the to- 
mogram is expected to vary as the rate of mass transfer from 
the normal companion into the Roche lobe of the white dwarf 
changes. 

The best agreement between the simulated and the ob- 
served tomograms was found assuming that the mass-transfer 
rate varies on a time scale of a few hours in the interval 
M ~ 10"' - 10'^ gs"' with the nightly mean value < M >^ 
5 X 10"' g s"' . In this case the efficiency of the propeller action 
by the white dwarf is limited to 0.01 ^ ist/^sd ~ 0.4. This 
means that the contribution of the propeller spin-down mech- 
anism to the observed braking of the white dwarf under the 
conditions of interest is small, and hence, the spin-down of the 
white dwarf (which is assumed to be governed by the pulsar- 
like spin-down mechanism) is expected to be stable indepen- 
dently of variations of M. 
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